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Abstract 

The precise driving force of the phase transition in indium nanowires on Si(lll) has been controversial 
whether it is driven by a Peierls instability or by a simple energy lowering due to a periodic lattice distortion. 
The present van der Waals (vdW) corrected hybrid density functional calculation predicts that the low- 
temperature 8x2 structure whose building blocks are indium hexagons is energetically favored over the 
room-temperature 4x1 structure. We show that the correction of self-interaction error and the inclusion of 
vdW interactions play crucial roles in describing the covalent bonding, band-gap opening, and energetics of 
hexagon structures. The results manifest that the formation of hexagons occurs by a simple energy lowering 
due to the lattice distortion, not by a charge density wave formation arising from Fermi surface nesting. 
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One-dimensional (ID) electronic systems have attracted much attention because of the richness 
of exotic physical phenomena such as charge density wave formation due to the Peierls instabil- 
ity II], non-Fermi liquid behavior |2], or Jahn-Teller distortion [3]. A prototypical example of 
quasi- ID systems is self-organized indium nanowires on the Si(l 1 1) surface |4J46|]. Each nanowire 
is composed of two zigzag chains of In atoms, and the nanowires are separated by a zigzag chain 
of Si atoms (see Fig. 1) [5]. Below ~120 K, this quasi- ID system undergoes a reversible phase 
transition initially from a 4 x 1 structure to a 4 x 2 structure, then to an 8 x 2 structure [ 6 , v|] , show- 
ing a period doubling both parallel and perpendicular to the In wires. This (4 x 1)^(8x2) phase 



transition is accompanied by a meta 



stimulated many experimental [|9l-ll2ll and theoretical studies 1113 



■insulator transition 



These intriguing results have 



42 ll| . However, the precise driv- 



ing force of the phase transition has been elusive for a long time. It has been suggested that the 
phase transition is driven b y a Peierls instability [|6l Isl — 1 1 ON or by a simple energy lowering due to a 



periodic lattice distortion [|13k 



2(m . The former mechanism involves the strong coupling between 



lattice vibrations and electrons near the Fermi level caused by Fermi surface nesting. Conse- 
quently, the charge density wave formation together with the lattice distortion occurs because of 



a larger electronic energy gain compared to an e 



mechanism involves either the trimer formation [113 



astic energy cost. On the other hand, the latter 



16fl in In chains with an elastic energy gain 



arising from the lattice distortion or the hexagon formation Ill7l - l2lll with an elastic energy gain 
from the lattice distortion as well as an electronic energy gain from the band-gap opening. 

(a) 




FIG. 1: Top view of the optimized (a) 4x1 and (b) 8x2 structures of the In/Si(lll) surface system. The 
dark and gray circles represent In and Si atoms, respectively. For distinction, Si atoms in the subsurface are 
drawn with small circles. Each unit cell is indicated by the solid line. 

Despite the above-mentioned controversial issue on the origin of the phase transition in the 
In/Si(l 11) system, the so-called hexagon model describes well several observed features of the 
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?D, scanning tunneling microscopy 



low-temperature phase such as an insulating character Q6J, |8 
images 111 in . and anisotropic optical interband transitions II 1211 . Initially, Gonzalez, Ortega, and 
Flores proposed that a shear distortion, whereby neighboring In chains are displaced in opposite 



1811 . Since this shear 



directions, allows for the formation of hexagon in the 4x2 unit cell IU7L 
phonon mode 111 ill is different from a phonon mode with the observed B6J] Fermi surface nesting 
vector 2kf = n/a x (a x is the lattice constant along the In chains), the Peierls mechanism is unlikely 
to be the driving force of the phase transition in the In/Si(l 1 1) system. Moreover, the stabilization 
of the 8x2 structure by doubling the unit cell perpendicular to the In wires is irrelevant with an 
electron-phonon coupling due to the observed Fermi surface nesting along the direction parallel 
to the In wires. The density-functional theory (DFT) calculations of Gonzalez, Ortega, and Flores 



showed that the 4x2 or 8 x2 hexagon structure is energetically favored over the 4x1 structure y/7 



2211 . but subsequent more accurate DFT calculations [1201 12 111 within the local density Il23ll as well 



as generalized gradient approximation VHAtt (LDA/GGA) did not support the energetic preference 
of the hexagon structures (see Table I). According to an LDA calculation with keeping the In 
Ad electrons frozen, the 8x2 hexagon structure was predicted to be energetically favored over 
the 4x1 structure D20J. However, this result is an artifact of the relatively less accurate scheme. 
Because of the proper prediction for the energy stability between the 4x1 and 8x2 structures, the 
LDA scheme with a frozen core of In Ad electrons has been employed to calculate the entropy 
difference [25] or the energy barrier between the two structures [26]. We note, however, that the 
LDA and GGA calculations with the treatment of the In Ad states as valence electrons predicted 
that the 4x2 and 8x2 hexagon structures are less stable than the 4x1 structure [|20j]. 

In this Letter, we present a new theoretical study which extends the previous work by consid- 
ering a hybrid exchange-correlation functional |[27ll and by taking van der Waals (vdW) Il28ll inter- 
actions into account. We will show that the correction of self-interaction error (SIE) cures over- 
delocalization of surface-state electrons inherent in the DFT and therefore describes adequately 
the covalent bonding, band-gap opening, and energetics of hexagon structures. Furthermore, we 
find that the vdW interactions between In atoms play an important role in further stabilizing the 
4x2 and 8x2 hexagon structures. Since the formation of hexagons and the more stabilization 
of the 8x2 structure are not associated with an electron-phonon coupling due to Fermi surface 
nesting, we can say that the phase transition in the In/Si(l 11) system is driven by a simple energy 
lowering due to the hexagon formation rather than by a Peierls-like mechanism. 

The present vdW corrected hybrid DFT calculations were performed using the FHI-aims I129ll 
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code for an accurate, all-electron description based on numeric atom-centered orbitals, with "tight" 
computational settings. For the exchange-correlation energy, we employed the hybrid functional 



of HSE [27] as well as the GGA functional of PBE OOfl . The k-space integrations in various 
unit-cell calculations were done equivalently with 64 k points in the surface Brillouin zone of the 
4x1 unit cell. The Si(l 1 1) substrate (with the Si lattice constant ao = 5.482 A) was modeled by 
a 6-layer slab (not including the Si surface chain bonded to the In chains) with ~15 A of vacuum 
in between the slabs. Each Si atom in the bottom layer was passivated by one H atom. All atoms 
except the bottom layer were allowed to relax along the calculated forces until all the residual 
force components were less than 0.02 eV/A. 

TABLE I: Calculated total energies (in meV per 4x 1 unit cell) of the 4x2 and 8x2 structures relative to the 



4x1 structure, together with the band gaps (in eV). For comparison, the previous LDA 02311 and GGA 02411 
results are also given. "Valence d (core d)" represents the treatment of the In 4d states as valence (core) 
electrons. 
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We begin to optimize the 4x1, 4x2, and 8x2 structures using the PBE functional. The op- 
timized 4x1 and 8x2 structures are displayed in Fig. 1(a) and 1(b), respectively. We find that 
the 4x2 and 8x2 structures show the formation of hexagons. Unlike previous pseudopotential 



calculations [13, 



14 



the present all-electron calculations were not able to find the stabiliza- 
tion of trimers in the 4x2 and 8x2 structures, which were converged to the 4x1 structure. The 
calculated total energies (in meV per 4x1 unit cell) of the 4x2 and 8x2 structures relative to the 



4x1 structure are given in Table I. We find that the 4x2 and 8x2 structures are less stable than the 
4x1 structure with A£4 X 2-4x1 = 33 meV and A.E8 X 2-4xi = 26 meV, respectively, consistent with 
those (48 and 27 meV in Table I) obtained by a previous GGA calculation ll20n with the Perdew- 
Wang exchange-correlation functional fl24|]. The calculated interatomic distances of In atoms are 
given in Table II. In the 4x1 structure, the In-In distance <ii ni -in 3 (^in 2 -in 4 ) within an In chain is 
3.045 (3.047) A, while din 3 -in 4 between the two In chains is 3.115 A. However, in the 4x2 (8x2) 
structure, the In-In distances between the two In chains are shortened as din 3 -in 4 = 3.054 (3.030) 
and Jin 7 -in 8 = 3.027 (3.011) A, leading to the formation of hexagon [see Fig. 1(b)]. We note that 
each In-In distance in the 8x2 structure slightly changes compared to the corresponding one in the 
4x2 structure because of the formation of hexagons in two opposite orientations (see Table II). 



TABLE II: Calculated interatomic distances (in A) of In atoms in the 4x1,4x2, and 8x2 structures using 
PBE. The results obtained using PBE+vdW are also given in parentheses. The labeling of In atoms are 
shown in Fig. 1. 
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To examine the influence of vdW interactions on the geometry and energetics, we use the 
PBE+vdW scheme developed by Tkatchenko and Scheffler Il28n. w here the vdW coefficients and 
radii are determined using the self-consistent electron density plQ . As shown in Table II, several 
In-In distances obtained using PBE and PBE+vdW exhibit some differences by less than 0.05 
A. The calculated PBE+vdW total energies of the 4x2 and 8x2 structures relative to the 4x1 



structure are also listed in Table I. We find that the 4x2 and 8x2 structures are still less stable than 
the 4x1 structure with A£4x2-4xi = 22 meV and A£gx2-4xi = 13 meV, respectively. Thus, the 
inclusion of vdW interactions within the PBE+vdW scheme does not reverse the stability of the 
4x2 (or 8 x2) and 4x1 structures. 

The calculated surface band structures of the 4x2 structure obtained using PBE and PBE+vdW 
are displayed in Fig. 2(a) and 2(b), respectively. It is seen that PBE gives a metallic feature 
while PBE+vdW gives an insulating feature with a band gap opening (Eg) of 0.05 eV. The PBE 
and PBE+vdW calculations for the 8x2 structure give almost the same value of E g = 0.08 eV 
[see Fig. 1(a) and 1(b) of the Supplemental Material (SI)]. For the 4x1 structure, both PBE and 
PBE+vdW predict well the observed metallic feature H, where three surface bands cross the 
Fermi level [see Fig. 2(a) and 2(b) of the SI]. Thus, we can say that PBE cannot predict the 
observed [6, 8, §| insulatin g fe ature for the low-temperature phase, consistent with a previous 



pseudopotential calculation 112 ill 
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FIG. 2: Surface band structure of the 4x2 structure obtained using (a) PBE, (b) PBE+vdW, (c) HSE, and 
(d) HSE+vdW. The inset in (a) shows the surface Brillouin zone for the 4x1 and 4x2 unit cells within that 
for the lxl unit cell. The energy zero represents the Fermi level. 



The local and semi-local DFT have challenged for a reasonable description of the energetics of 
different structures involved in the phase transition or the kinetics of chemical reactions because 
of their intrinsic SIE [32, 2^. 



Especially, the GGA tends to stabilize artificially delocalized states 



due to the SIE, since delocalization reduces the self-repulsion. It is thus likely that the present 
PBE functional would give a lower energy for the metallic 4x1 structure, compared to the 4x2 
and 8x2 structures. In order to correct the SIE, we use the hybrid HSE functional [|27j]_to calculate 
the total energies of the 4x1, 4x2, and 8x2 structures with the PBE geometries (13411 . We find 
that the correction of SIE enhances the stability of the 4x2 and 8x2 structures relative to the 



4x 1 structure, giving rise to A£4 X 2-4xt = 3 meV and A£8x2-4xi = — 15 meV, respectively. This 
enhanced HSE stability of the 4x2 (8x2) structure is caused by the electronic energy gain arising 
from an increased band gap of E g = 0.10 (0.19) eV, as shown in Fig. 2(c) [Fig. 1(c) of the SI]. 
Thus, HSE predicts well the observed DalSllSD insulating feature for the 4x2 and 8x2 structures. 

To see the effects of the SIE on the charge density distribution, we plot the charge density 
difference defined as 

Ap = Phse-Ppbe, (1) 

where Phse (Ppbe) is the charge density obtained using the HSE (PBE) functional. The results for 
the 4x1 and 4x2 structures are displayed in Fig. 3(a) and 3(b), respectively. We find a conspicuous 
difference between the 4x1 and 4x2 structures for Ap. It is seen that the insulating 4x2 structure 
has a larger Ap between In atoms compared with the metallic 4x1 structure, indicating that the 
relatively localized surface states in the former are more affected by the SIE than the delocalized 
surface states in the latter. This fact also reflects that in the 4x2 structure, the correction of the 
SIE by HSE recovers the charge localization in the covalently bonding between In atoms. As 
a consequence, we obtain an increase of band gap with E g = 0.10 eV, leading to a decrease of 
A£4x2-4x l = 3 meV compared to the PBE result (A£ , 4 X 2-4 X l = 33 meV). For the 8x2 structure, 
Ap shows a similar pattern with the 4x2 case (see Fig. 3 of the SI), yielding E g = 0.19 eV and 
A£8x2-4x l = — 15 meV. We note that there is a general trend that the 8 x2 structure is more stable 
than the 4x2 structure (see Table I). This indicates some energy gain caused by the correlation 
between two In nanowires in the 8 x 2 structure, as pointed out by a previous theoretical study [20] . 



(a) 



(b) 





FIG. 3: (Color on line) Charge density difference between Phse and Ppbe for the (a) 4xf and (b) 4x2 
structures. The dark (gray) color represents the isosurface of 0.02 (—0.02) e/A 3 . 



Using the HSE+vdW schme, we ca 
tures with the PBE+vdW geometries 



culate the total energies of the 4x1,4x2, and 8x2 struc- 
3411 . We find that the 4x2 and 8x2 structures are more 
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stable than the 4x1 structure with A£4 X 2-4xi = —23 meV and A.E8x2-4xi = —40 meV, respec- 
tively. Since the total energy is composed of the HSE energy (£hse) and the vdW energy (£Vdw) 
which is given by a sum of pairwise interatomic C(,R 6 terms, the total energy difference between 
the 4x2 (or 8x2) and 4x1 structures is determined by 



A£ = A£ H se + A£ 



WdW- 



(2) 



Figure 4 shows AfiVdw m AE4x2-4xi and AE^ X 2-4xi, together with its components originating 
from In— In, In— Si, and Si— Si atoms. We find that A£hse in AE4 X 2-4xl (A£8x2-4xi) is —4 (—18) 
meV, while A.E v dw in AE4 X 2-4x l (A£8x2-4x l) is — 19 (—22) meV. Therefore, the inclusion of vdW 
interactions largely enhances the stabilization of the 4x2 and 8x2 structures. We note that the 
A£hse values (—4 and —18 meV) in A£4 X 2-4x i and A£8 X 2-4x l are somewhat different from those 
(3 and —15 meV) obtained from the HSE calculation due to the use of two different PBE+vdW 
and PBE geometries. As shown in Fig. 4, the component of AZsVdw originating from In— In atoms 
in the 4x2 and 8x2 structures is significantly larger in magnitude than those originating from 
In— Si and Si— Si atoms. Thus, we can say that the vdW interactions between In atoms play an 
important role in stabilizing the formation of hexagons. 
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FIG. 4: (Color on line) Contribution of vdW energy to the total-energy difference A£4 X 2-4 X i (AZ?8 X 2-4xi) 
obtained using HSE+vdW. The components originating from In— In, In— Si, and Si— Si atoms are also given. 



As shown in Fig. 2(d) [Fig. 1(d) in the SI], the HSE+vdW band structure of the 4x2 (8x2) 
structure gives E g = 0.21 (0.21) eV, in good agreement with a recent scanning tunneling spec- 
troscopy measurement of E g = 0.2 eV (13511 . We note that the HSE+vdW magnitude of A.E8x2-4x l 
is 40 meV per 4x1 unit cell, which is equal to 10 meV per In atom. Although the precise 
(4x l)<f+(8x2) phase transition temperature can be estimated by comparing the vibrational free en- 



ergies of the 4x1 and 8x2 structures [25], the HSE+vdW magnitude of A£8 X 2-4x l is well compa- 
rable with the thermal energy at the observed phase transition temperature of ~120 K | (], ^]. There- 
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fore, the HSE+vdW scheme is likely to give good band gap and energetics of the low-temperature 
phase in the In/Si(l 1 1) system. 

In summary, using the HSE and HSE+vdW schemes, we investigated the energy stability of the 
low-temperature and room-temperature structures in the In/Si(lll) system, which has not been 



2 ill . We found that the correction of SIE 



adequately described by previous DFT calculations H20L 
cures the derealization error not only to give the insulating feature for the 4x2 and 8x2 structures 
but also to reverse the stability of the 4x1 and 8x2 structures. We also found that the vdW inter- 
actions between In atoms enhance the stability of hexagon structures. Our results demonstrate that 
the formation of hexagons in the In/Si(l 11) system occurs by a simple energy lowering due to the 
lattice distortion rather than by a Peierls instability. We notice that the Sn/Si(l 1 1) and Sn/Ge(l 1 1) 
surfaces have been the object of a large number of studies for determining the exact crystallo- 
graphic arrangement, the electronic structure, and the mechanism of the phase transition I136tl . 
We anticipate that the correction of self-interaction error and the inclusion of vdW interactions 
would give more accurate description for the structural and electronic properties of such prototype 
two-dimensional electron systems. 

This work was supported by National Research Foundation of Korea (NRF) grant funded by 
the Korean Government (NRF-201 1-0015754). The calculations were performed by KISTI super- 
computing center through the strategic support program (KSC-2012-C3-11) for the supercomput- 
ing application research. 
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